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Abstract 

By exhibiting the corresponding Lax pair representations we propose a wide class of integrable 
two-dimensional (2D) fermionic Toda lattice (TL) hierarchies which includes the 2D N = (2|2) 
and N = (0|2) supersymmetric TL hierarchies as particular cases. Performing their reduction 
l— ~~ to the one-dimensional case by imposing suitable constraints we derive the corresponding ID 
y—i ■ fermionic TL hierarchies. We develop the generalized graded R-matrix formalism using the 
generalized graded bracket on the space of graded operators with an involution generalizing 
■^j- . the graded commutator in superalgebras, which allows one to describe these hierarchies in the 
framework of the Hamiltonian formalism and construct their first two Hamiltonian structures. 
The first Hamiltonian structure is obtained for both bosonic and fermionic Lax operators while 
the second Hamiltonian structure is established for bosonic Lax operators only. We propose the 
graded modified Yang-Baxter equation in the operator form and demonstrate that for the class 
of graded antisymmetric i?-matrices it is equivalent to the tensor form of the graded classical 
\ Yang-Baxter equation. 
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1 Introduction 



The Toda lattice (TL) is one of the most important families of models in the theory of integrable 
systems. Its various generalizations and supersymmetric extensions, having deep implications 
in modern mathematical physics, have been the subject of intense investigations during the last 
decades. 

The 2D TL hierarchy was first studied in PJ Ej, and at present two different nontrivial 
supersymmetric extensions of 2D TL are known. They are the N = (2|2) 0-^3 and N = (0|2) 
[T2"l |§| supersymmetric TL hierarchies that possess a different number of supersymmetries and 
contain the N = (2|2) and iV = (0|2) TL equations as subsystems. Quite recently, the 2D 
generalized fermionic TL equations have been introduced ^T] and their two reductions related 
to the N = (2 1 2) and N = (0|2) supersymmetric TL equations were considered. In the present 
paper, we describe a wide class of integrable two-dimensional fermionic Toda lattice hierarchies 
- 2D fermionic (K + ,K-)-TL hierarchies, which includes the 2D N = (2|2) and N = (0|2) 
supersymmetric TL hierarchies as particular cases and contains the 2D generalized fermionic 
TL equations as a subsystem. 

The Hamiltonian description of the 2D TL hierarchy has been constructed only quite re- 
cently in the framework of the R-matrix approach in [T^J, where the new R-matrix associated 
with splitting of algebra given by a pair of difference operators was introduced. In the present 
paper, we adapt this R-matrix to the case of ^-graded operators and derive the bi-Hamiltonian 
structure of the 2D fermionic (K + , K~)-TL hierarchy. 

Remarkably, in solving this problem the generalized graded bracket (j5J) on the space of 
graded operators with an involution finds its new application. This bracket was introduced in 
[T3] . where it was observed that the N — (1|1) supersymmetric 2D TL hierarchy had a natural 
Lax-pair representation in terms of this bracket which allowed one to derive the dispersionless 
iV = (1|1) 2D TL hierarchy and its Lax representation. In the present paper, the generalized 
graded bracket is used to describe the 2D fermionic (K + ,K~)-TL hierarchy and define its 
two Hamiltonian structures. Moreover, we demonstrate that the classical graded Yang-Baxter 
equation J3| has an equivalent operator representation in terms of the generalized graded 
bracket. All these facts attest that this bracket has a fundamental meaning and allows a broad 
spectrum of applications in modern mathematical physics. 

This paper is the extended version of [IE] . The structure of the paper is as follows. 

In Sec. 2, we define the space of the Z 2 -graded difference operators with the involution and 
recall the generalized graded bracket ^3] and its properties. 

In Sec. 3, we give a theoretical background of the .R-matrix method generalized to the case 
of the ^-graded difference operators. We define the .R-matrix on the associative algebra g of 
the Z 2 -graded difference operators, derive the graded modified Yang-Baxter equation and using 
the generalized graded bracket obtain two Poisson brackets for the functionals on gt = g. The 
proper properties of the Poisson brackets thus obtained are provided by the properties of the 
generalized graded bracket. Thus, for the ^-graded difference operators of odd (even) parity 
this bracket defines odd (even) first Poisson bracket. The second Poisson bracket is found only 
for even difference operators which in this case are compatible with the first Poisson bracket. 
Using these Poisson brackets one can define the Hamiltonian equations that can equivalently 
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be rewritten in terms of the Lax-pair representation. The basic results of Sec. 3 are formulated 
as Theorem. 

In Sec. 4, using the generalized graded bracket we propose a new 2D fermionic (K + ,K~)- 
TL hierarchy in terms of the Lax-pair representation and construct the algebra of its flows. 
Then we present the explicit expression for its flows and show that all known up to now 2D TL 
equations can be derived from this hierarchy as subsystems. 

In Sec. 5, we consider the reduction of the 2D fermionic (K + ,K~)-TL hierarchy to the 
ID space and reproduce the ID generalized fermionic TL equations [llj as the first flow of the 
reduced hierarchy with additional constraint imposed. 

In Sees. 6 and 7, we apply the results of Sec. 3 to derive the Hamiltonian structures of the 
ID and 2D fermionic (K + , i^ _ )-TL hierarchies. Following ^3] we use the i?-matrix which acts 
nontrivially on the space of the direct sum of two difference operators and derive two different 
Hamiltonian structures of the 2D fermionic (K + ,K~)-TL hierarchy. The first Hamiltonian 
structure is obtained for both even and odd values of (K + , K~) while the second one is found 
for even values of (K + ,K~) only. We perform their Dirac reduction and demonstrate that in 
general the Dirac brackets for the second Hamiltonian structure are nonlocal but for the case of 
the fermionic (2, 2)-TL hierarchy they become local. As an example, we give the explicit form 
of the first and second Hamiltonian structures for the fermionic ID (2,2)-TL hierarchy. 

In Sec. 8, we briefly summarize the main results obtained in this paper and point out open 
problems. In Appendices, we clarify some technical aspects. 



2 Space of difference operators 

In this section, we define the space of difference operators which will play an important role in 
our consideration. These operators can be represented in the following general form: 

oo 

°™ = E /i!?^ 9 , mjez, (1) 

k=— oo 

parameterized by the functions f^j (fit+ij) which are the Z 2 -graded bosonic (fermionic) lattice 
fields with the lattice index j (j e Z) and the Grassmann parity defined by index k 

dj m ) = \k\ mod 2. 

Jk,j 

e kd is the shift operator whose action on the lattice fields results into a discrete shift of a lattice 
index 

Jd Am) _ Am) Id / 9 \ 
e Jk,j — JkJ+l^ ■ \ z ) 

The shift operator has Z 2 -parity defined as 

d'w = \l\ mod 2. 
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The operators O m (JTJ) admit the diagonal Z 2 -parity 



and the involution 



d® m = d jm) + d' e(k _ m)a = \m\ mod 2 (3) 

■* k .•' 



E (- 1 W e( 



^fe f{m) e (k-m)d 

"m 

k=— oo 



In what follows we also need the projections of the operators O m defined as 



ki^p+m k^p+m 



e (k-m)d 



and we will use the usual notation for the projections (Q m )+ := (O m )^ an d (O m )_ '■= 
(O m )<o- Note that e ld is a conventional form for the shift operators defined in terms of infinite- 
dimensional matrices (e ld )ij = dij-i, and there is an isomorphism between operators (JIJ and 
infinite-dimensional matrices (see e.g. [T7j ) 

oo oo 

®m= E /£? e(fc " m)9 - - E /£? 

fc=— oo fc=— oo 

In the operator space (JTJ) one can extract two subspaces which are of great importance in 
our further consideration 



fc=0 



E^ ±(K± "* )8 . ^ ±eN - ( 4 ) 



The operators of the subspaces form associative algebras with the multiplication (J2J). Using 
this fact we define on these subspaces the generalized graded algebra with the bracket 

[O, 6} :=00- (-l) d o d o Q*H>) 0* (d 6 } , (5) 

where the operators O and O belong to the subspaces 0^ + ((QC_), and 0**" 1 ) denotes the 
m-fold action of the involution * on the operator O, (O*^ 2 - 1 = O). Bracket (jSJ) generalizes the 
(anti) commutator in superalgebras and satisfies the following properties 
symmetry 

[O, O} = -(-l)<Vd [Q*( d o) ; 0* (d o)}, (6) 

derivation 

[o,od} = [o,6}6 + (-i)^6*W[o*w,6}, (7) 

and Jacobi identity 

+ (_i)Vd [[6, 0* (d 6 } }, d* (d & +cto) } = 0. (8) 
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For the operators O m (JTJ) we define the supertrace 

oo 

strO= ^(-l)VS- (9) 

j=-oo 

In what follows we assume suitable boundary conditions for the functions f& in order the 
main property of supertraces 

str[O,O} = (10) 
be satisfied for the case of the generalized graded bracket (J5J). 



3 R-matrix formalism 

In this section, we develop a theoretical background of the R-matrix method adapted to the 
case of the operator space (JTJ) . 

Let g be an associative algebra of the operators from the space (P) with the invariant 
non-degenerate inner product 

< 0,6 >= str (OO) 

using which one can identify the algebra g with its dual g^. We set the following Poisson 
bracket: 

{/,<?}(©) = - < O, \Vg, (V/)*^} >, (11) 

where /, g are functionals on g, and V/ and Vg are their gradients at the point O which are 
related with /, g through the inner product 



df(0 + e 



de 



< SO, V/(0) > . 



e=0 



Note that the proper properties of the Poisson bracket (|11|) follow from the properties (J0HSI) 
of the generalized bracket (0) and are strictly determined by the Z 2 -parity of the operator O. 
Thus, one has symmetry 

{f,g} = -(-l)(df+^ B +^{gJ}, (12) 

derivation 

{f,9h} = {f,g}h+(-l) d ^ d f + ^g{f,h}, (13) 

and Jacobi identity 

(_ 1 )(rf/+*>)(«k-Hio){{£ g ^ h y + (_ 1 )K+do)(rf/+rfo){{^ ) h ^ 
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Therefore, for the even operator O one has a usual (even) ^-graded Poisson bracket, while 
for the operators with odd diagonal parity d Q (J2J) eq. (fTTj) defines the odd Z 2 -graded Poisson 
bracket (antibracket). 

Having defined the Poisson bracket we proceed with the search for the hierarchy of flows 
generated by this bracket using Hamiltonians. Therefore, we need to determine an infinite 
set of functionals which should be in involution to play the role of Hamiltonians. For Poisson 
bracket (jllj) one can find an infinite set of Hamiltonians in a rather standard way 

oo 

H k = ^strOt = ±J2 (-!)'/£#, ( 15 ) 

i=— oo 

where is defined as 

(Q)f := (0* (do) 0) fc , (O)f +1 := O (O)f . (16) 

For the odd operators O eq. (J15)) defines only fermionic nonzero functionals i^fe+i, since in 
this case even powers of the operators O have the following representation: 

d = 1 : (O)f = (l/2[(0)*, 0}) k = l/2[((0)f - 1 )*, O} (17) 

and all the bosonic Hamiltonians are trivial (i?2fc = ) like the supertrace of the generalized 
graded bracket. 

The functionals 1)15)1 are obviously in involution but produce a trivial dynamics. Actually, 
the functionals (|15j) are the Casimir operators of the Poisson bracket (|11|). so the Poisson 
bracket of with any other functional is equal to zero as an output (due to the relation 
Viffe+i = O*). Nevertheless, it is possible to modify the Poisson bracket ((TT} in such a way 
that the new Poisson bracket would produce nontrivial equations of motion using the same 
Hamiltonians (fTK)) and these Hamiltonians are in involution with respect to the modified Poisson 
bracket as well. Let us introduce the modified generalized graded bracket on the space (JTJ 

[0,€>} R := [R{0),®} + [0,R{0)}, (18) 

where the .R-matrix is a linear map R. g — > g such that the bracket (|T%j) satisfies the properties 
(JSHHl- One can verify that the Jacobi identities (jHJ) for the bracket (JTHj) can equivalently be 
rewritten in terms of the generalized graded bracket (j3J) 

[[o, 0*^} fl , d< d o+ d o)j R + cyc i e perm. = 
[^([o, 6* {d o)} fl ) - [i?(0), ^(O*^)},©*^^} + cycle perm. = 0. 

Thus, one can conclude that a sufficient condition for R to be the i?-matrix is the validity of 
the following equation: 

R([D, 0} R ) - [R(0), R(D)} = a[Q, 6}, (19) 

where a is an arbitrary constant. One can show (see Appendix A) that for the case of graded 
antisymmetric operators R eq. ()19)1 at a = 1 represents the operator form of the graded 
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classical Yang-Baxter equation ^H]- Following the terminology of ^H] we call eq. (fTTI|) the 
graded modified Yang-Baxter equation. Equation (fTUj) is the generalization of the graded 
modified classical Yang-Baxter equation discussed in paper [TH] for the case of space of graded 
operators (JTJ). 

With the new bracket ()18|) one can define the corresponding new Poisson bracket on dual g^: 
{f,g}i(0) = - 1/2 < O, [V<?, (V/)*^ fl )} B > 

= I < (_l) rf V 9 do^yg)*(do)p*(dvg) ; (\/fy( d Vg) j 

- [0,Vg}R((Vf)< d ^)>. (20) 

With respect to the dependence of the r.h.s of (j2*Uj) on the point O this is a linear bracket. 
Without going into details we introduce also bi-linear bracket for bosonic graded operators O b 
( do B — 0) as follows: 

{/,^} 2 (Ob) = -l/4<[0 B ,Vg}R((Vf)< d ^0f^ +d ^ 
+ 0*^ s) (Vf)* id ^) - R(\7gQ^ dv9) 

+ B Vg)[0^ dv3 \ (V/)*^} > . (21) 

We did not succeed in constructing the bi-linear bracket for the case of fermionic operators 
( do F = 1). The bracket (j2*U|) is obviously the Poisson bracket if R is an i?-matrix on g. The 
bi-linear bracket (j21j) becomes Poisson bracket under more rigorous constraints which can be 
found in the following 

Theorem, a) Linear bracket (|2T)|) is the Poisson bracket if R obeys the graded modified 
Yang-Baxter equation fTT)j) : 

b) the bi-linear bracket (|21j) is the Poisson bracket if R and its graded antisymmetric part 
l/2(i? — R^) obey the graded modified Yang-Baxter equation (JTI5J) with the same a; 

c) if O = Ob, then these two Poisson brackets are compatible; 

d) the Casimir operators Hk |T5|) of the bracket (|TT| are in involution with respect to both 
linear (|2T)|) and bi-linear (|2*T]) Poisson brackets; 

e) the Hamiltonians Hk ^ (|TK|) generate evolution equations 

d k O = {i? fc+1 ,0} 1 = l/2[ J R((ViJ fc+ i)* (do) ) ) 0}, 
d k O B = {H k ,Q B } 2 = l/4[R(VH k Q B + B VH k ),0 B } (22) 

via the brackets (J2(Jj) and (J21|) . respectively, which connect the Lax-pair and Hamiltonian rep- 
resentations. 

Proof, a). This is a summary of the above discussion on the linear bracket. 

b) . Using the property of cyclic permutations inside the supertrace one can easily verify that the 
symmetry property {/, g} 2 = —(—{) dfd9 {g,f}2 holds. Verification of the Jacobi identities for 
the bi-linear bracket amounts to straightforward and tedious calculations which are presented 
in Appendix B. 

c) . These two Poisson brackets are obviously compatible. Indeed, a deformation of the point 
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Ob — > Ob + b on the dual g^, where b is an arbitrary constant operator, transforms (j21|) into 
the sum of two Poisson brackets 

{/, 9h(®B + b) = {/, gh(o B ) + b{f, gUO B ). 

d) . Substituting the expressions for Hamiltonians (fTKj) into eqs. (ffilj) and (J2IJ and taking into 
account that V-fffc+i = it is easily to check that the Casimirs of the bracket are in 
involution with respect to both the Poisson structures (|2"U|) and (}2"T]) . 

e) . Using cyclic permutations inside the supertrace operation let us rewrite both the Poisson 
brackets (|20|) and (|21|) in the following general form: 

{/, gUO) =< Pi(p) Vg, (V/)*^ >, i = l, 2, 

where P{(0) is the Poisson tensor corresponding to the bracket {..,..}, 

P 1 (0)Vg = -l/2([0,R(Vg)} + R\[0,Vg}})), 
P 2 (0 B )Vg = l/A([R(VgOB + B Vg),0 B {dv3) } 

- B R\[0 B , Vg}) - R\[O b , Vg})Of V9) ) 

and the adjoint operator RJ acts on the dual g^ 

< Q,R(b) >=< R\0),6> . 

The Hamiltonian vector field associated with Hamiltonian is given by dkO = Pj(0)Viffc. 
Taking into account that [O, VH^} = we arrive at the Lax-pair representations (|22|) . ■ 

Note that a similar Theorem when the shift operators and functions parameterizing the 
difference operators O (JTJ) have even Z 2 -parity was discussed in [THJ I20J OH HI] ■ 

For the graded modified Yang-Baxter equation (j!9|) there is a particular class of solutions 
which are useful in application. Suppose that the algebra g can be represented as a vector space 
direct sum of two subalgebras 

g = g+ + g-: [g +! g+}Cg + , [g_,g_}cg_. 

Let P± be the projection operators on these subalgebras, P±g = g±, then one can easily verify 
that R = P + — P- satisfies the graded modified Yang-Baxter equation (fT^j) at a = 1 and, 
therefore, represents the i?-matrix on g. Indeed, in this case the modified generalized graded 
bracket (HHJ) 

[0, 0} R = 2[(0)+, (6)+} - 2[(0)_, (6)_} (23) 

obviously satisfies the Jacobi identities (jHJ), since it determines the usual direct sum of two 
subalgebras 

[g±,g±}i? C g±, [g + ,g-} R = 0. 
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4 2D fermionic (K + 1 K )-Toda lattice hierarchy 

In this section, we introduce the two-dimensional fermionic (K + , K~)-Toda lattice hierarchy in 
terms of the Lax-pair representation. 

Let us consider two difference operators L K ± 

OO OO 

L + K+ = £ W^-W L~_ = £ v^-V, (24) 

k=0 k=0 

which obviously belong to the spaces (@J). The lattice fields and the shift operator entering into 
these operators have the following length dimensions: [u^i] = —l/2k, [vk t i] = l/2(k — K + — K~) 
and [e kd ] = —l/2k, respectively, so operators are °f equal length dimension, [L~L + ] = 
[L"Z_] = —1/2K + . The dynamics of the fields u^^v^^ are governed by the Lax equations 
expressed in terms of the generalized graded bracket (0) jE] 

DfL% a = T a(-ir KaK± mL^ ± )l)^)< Ka \L a Ka }, <* = +,-, s G N, (25) 

where Df are evolution derivatives with the ^-parity defined as 

d D ± = sK^ mod 2 

and the length dimension [Df] = [D~] = —sK + /2. The Lax equations ()25|) generate non- 
Abelian (super) algebra of flows of the 2D fermionic (K + , K~)-TL hierarchy 

[Df, D±} = (1 - (-ir K± )Df +p , [D+, D;} = 0. 

The composite operators (i/5. ± )* entering into the Lax equations (|2*3jl are defined by eq. (fTfij) 
and also belong to the spaces (jlj) 

OO OO 

k=0 k=0 

Here and v£\ are functionals of the original fields and there are the following recursion 
relations for them 

v 

(r+l) _ \^(_ ]) kK+ (r) 7 ,W-„ • 

k=0 

V 

(r+l) _ (r) 7, (1) -7, . 

%V ~~ ' V k,i V p-k,i+k-rK~, V p>i —V Pll . 

k=0 

Now using the Lax representation (|2*5|l and relations (J7|) and (fT^j) one can derive the equations 
of motion for the composite Lax operators 

Df{Lu: = ^(-i) s ^^ ± [(((L| ± ):)-a)* m ,(^);}. (26) 
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The Lax-pair representation (|2T)|) generates the following equations for the functionals 



(r) (r) 



Dt4} = j2((- l ) rpK++1JS) "" j 



p+sK+,i k— p,i— p 
p=\ 



i (_l)(k+p)sK+ (r) (a) \ ^ 

\ L ) a k- v ,i a r ) +sK+.i+v-k+rK+h \*') 



L k-p,i Uj p+sK+,i+p-k+rK+J' 
sK--l 

D -JX) _ ( (_ 1 \(sK-+p)rK+ (a) (r) 

U s U k,i — 2-^1 ^ > U p,i U p+k-sK-,i+p-sK- 
p=0 



,s(k+p+l)sK- (r) (s) \ fns 

-V a v+k-sK-.i u v.i-v-k+sK-+rK+ h 



p+k—sK~,i p,i—p—k+sK~+rK+'i 
sK+ 

p=0 

_ (_l)(k+P+l)sK+ (r) (*) v , 2q ^ 

V x / v p+k-sK+,i a p,i+p+k-sK+-rK~ ■>> V^ 3 7 

i C_lUfc+p)5-K- (r) (a) n / oqN 

+ ^ ^ U k-p,i U p+sK-,i+k~p-rK-)- 

It is assumed that in the right-hand side of eqs. (|27H3U|l all the functionals v/[\ with k < 
should be set equal to zero. 

Let us demonstrate that all known up to now 2D supersymmetric Toda lattice equations 
can be derived from the system (j2?H30j) . 

First, the 2D generalized fermionic Toda lattice equation discussed in [TT] can be reproduced 
from the system of equations (|2*THHUj) as a subsystem with additional reduction constraints 
imposed. In order to see this, let us introduce the notation t> 0) j = di, v\ t i = pi, u\ t i = 7$, «2,« = Q 
and consider eqs. (|2*5|) and (|50|) at K + = K~ = 2, r = s = 1. One obtains 



Dfd; = di(ci - q_ 2 ), £>i~7i = PiU 0ii -i - Pi+2«o,i 

Cj = diU 0ti _ 2 ~ d i+2 uo t i - liPi+i - ji-iPi, 
Dfpi = pi(ci- + di+i7i- diji-2, D^u 0ti = 0. (31) 

It is easy to check that after reduction uo,i = 1 eqs. ffTTTf) coincide with the 2D generalized 
fermionic Toda lattice equations up to time redefinition D± — > —D^. 

Next, the N=(2|2) supersymmetric Toda lattice equation also belongs to the system (|2*7H5n|) . 
In order to see that, let us consider eqs. ()28|) at K + = K~ = s = k = r = l 

D^uij = -u ,i-iv ,i - u ,iV 0ii+ i (32) 
and eqs. (12H at K + = K~ = s = r = 1. k = 

£*l"Ho,< = Uo,i(wi,i - (33) 
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Then imposing the constraint uo,i = 1 and eliminating the fields from eqs. ()H2fJ3*3*|) one 
obtains the N=(l|l) superfield form of the N=(2|2) supersymmetric Toda lattice equation 

0001121 

D\D^ In v ,i = f o,i+i - v ,i-i. (34) 

Analogously, one can show that the consideration of eqs. ()29j) and (j50"|) at K + = 1, K~ = 2, 
s = r = 1 and k = 0,1 leads to the N=(0|1) superfield form of the N=(0|2) supersymmetric 
Toda lattice equation [121 E] after imposing the reduction constrains uo,« = 1, fo,2i+i = 0. 

We call equations ([251) f° r arbitrary (K + ,K~) the 2D fermionic (K + , K~)-Toda. lattice 
hierarchy. 



5 ID fermionic K )-Toda lattice hierarchy 

In this section, we consider the reduction of the 2D fermionic (K + , K")-Toda lattice hierarchy 
for even values of (K + , K~) to the ID space. 

Let (K + ,K~) be even numbers. In this case the generalized graded bracket (0) between 
two Z2-even operators turns into the usual commutator and eqs. ()25|) become 

DfL a Ka = [((L± ± y) ± ,LU- (35) 

Following [22] for even (K + , K~) one can impose the reduction constraint on the Lax oper- 
ators (J24|) as follows: 

L + K+ + (L+ + )" 1 = L~_ + (L-_)- x = L K+tK -, (36) 
which leads to the following explicit form for the reduced Lax operator 

K+ K~-l K++K- 

L K+ , K - = J2^ {K+ ' k)d + E v^ k - K ~ )9 = E ^e [K+ - k)d - (37) 

k=0 k=0 k=0 

Substituting the expressions for the reduced composite Lax operators L K± = L K + K - — (L^ ± ) _1 
into Lax equations (|33|) one can see that these equations become equivalent to the single Lax 
equation on the reduced Lax operator 

D S L K + )K - = l((L K + )K -) s ) + ,L K + yK -} (38) 

with = —D~ — D s . As a consequence of eq. (f3~%|) we have 

D s (L K + tK -) r = [((L K + jK -) s ) + ,(L K + iK -) r ). 

At the reduction uo,i = 1 the ID (2, 2)-TL hierarchy becomes that studied in detail in [TT] . 
In this case, the representation (j3*H)) with the Lax operator 

L 2 , 2 = e 29 + 7i e 9 + a + pte' 9 + die' 29 
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gives the following first flow: 

D\di = di(ci — Q-2)) 

DxPi = Pi{ci - Cj_i) + di+x'ji - diji-2, 

Dlji = Pi+2 - Pi, 

D x Ci = d i+2 -di + 'jiPi+i + li-iPi- 

These are the ID generalized fermionic Toda lattice equations ^T] which possess the N = A 
supersymmetry. 

6 Bi-Hamiltonian structure of ID fermionic i^~)-TL 
hierarchy 

In this section, we apply the R-matrix approach to build the bi-Hamiltonian structure of the 
ID fermionic (K + , K~)-TL hierarchy and perform its Dirac reduction. 

The space of operators Q^.+ (HJ) can obviously be split into the vector space direct sum, 
Ok+ = (Ok+)+ + (Ok+)— The .R-matrix arising from this splitting 



R = P + -P_, i?(Oi 



»K+ + ~ 



obviously solves the graded modified Yang-Baxter equation (j!9|) at a = 1. This R- matrix is not 
graded antisymmetric, R 7^ — R?\ however, its graded antisymmetric part A = 1/2(R — R^) as 
well as the .R-matrix itself satisfy the graded modified Yang-Baxter equation (J19|) . According to 
the general Theorem of Section 3 this means that there exist two Poisson structures on gt = g. 
Substituting the general form of operators L~^ + (|2"2j) and 



Vu n4 = e^- K+ ^ d (-iy5^ 

into (|2(J|) and (J21j) one can find the explicit form of the first and second Poisson brackets, 
respectively, 

{u n ,i,u m j}i = (-iy(5~ K+ +5 mK+ -l)(u n+rn _ K+ji 5 iij+n _ K + 

I 1 \(m+K + )(n+K+ + l) n r / Q rA 



and 



n+m 



k=0 



(_l)m(rH-fc+l) Ufc AUn+m _ k j S i; j _ m+k )] , (40) 
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where 

+ f 1, if n > m - ( 1, if n < m 

n ' m = \ 0, if n < m ' n > m = \ 0, if n > m. 

Let us remind that the second Poisson brackets are defined for even values of K + only. 

Our next goal is to perform the reduction of Poisson brackets (|39H4(Jj) for the functions 
parameterizing the operators L~j <+ (J24j) to the Poisson brackets corresponding to the reduced 
operators (JHTj) 

k++k- 

k=l 

where K + , K~ are even numbers. Therefore, one needs to modify the Poisson brackets 
according to the reduction constraints 



Uk i = 0, k > K + + K~ 
u ,i = 1 



for any i. (41) 



We apply these reduction constraints in two steps. First, we note that for the first constraint 
in ()41|) the reduction simply amounts to imposing constraint w^j = (k > K + + K~) due to 
the observation that 



« M =o, =0, 0<n<K + + K-, m>K + + K-, p = 1,2. (42) 

k>K++K~ 



For the first Poisson brackets (JHUj) relation (J4*2")) is obvious. One can derive eq. (J4*2*|) for the 
second Poisson brackets (j4*Uj) if one divides the sum in (J4*U|) into three pieces 

n + m max(0.n+m-K+-ft''-l) vaai(n+m-l,K + +K~) n+m 

E = E + E + E («) 

k=0 k=0 k=ma,x(l,n+m-K+-K-) k=mm(n+m,K++K-+l) 

Now it is easy to verify that the second sum in the r.h.s. of eq. (j4*Hj) is the only sum which 
could give a nonzero contribution to eq. (j4lZ|) , but it is equal to zero if < n < K + + K~ , m > 
K+ + K-. 

Now let us consider the second reduction constraint in eq. (j4*Tj) . w 0) i = 1. Following the 
standard Dirac reduction prescription we obtain for the Dirac brackets 



«0,i=l 



with the correction term 



i ,3 



U0,i=l 
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In the case of the first Poisson brackets {uo,i, u m j}\ = for any m. Thus, one can conclude 
that the first Poisson brackets (J3TJJ) are not modified, and they can simply be restricted by 
imposing the constraints (|4ip. 

Before investigating the second Poisson brackets (jjUj) we supply the fields u n ^ and v n ^ with 
the boundary conditions 

lim u U) i = lim v n ^ = 0, (44) 

and introduce a new notation 5ij- n := {A n )ij which is useful in what follows. One can verify 
that in the new notation the multiplication of matrices results in adding powers of the operators 
A: 

Then, one can represent the correction term for the reduced second Poisson brackets as follows: 

A 2 («n,i, Um,j) = l/2(-l)V ! [(l-A Jf+ -')(l + A- K+ ) 

(A* + - A~ K+ )-\l + A* + )(l - (-l) m A m - K+ )] hJ u mj . 

In general, the reduced second Poisson brackets are nonlocal since the inverse matrix 

(A 2u - A' 2 ")- 1 = (1 + A 2v )-\l - A- p )-\l + A"")- 1 , 

being considered as an operator acting in the space of functionals with boundary conditions 
can be expressed via infinite sums 



oo oo 

'l-hr v )- x = A x ^A- fel/ - (1 - AOj^A*", 

k=0 k=l 

oo oo 

[1 + A-")- 1 = A 2 ^(-l) fc A-^ + (l-A 2 )^(-l) fe A^, 

fc=0 fc=l 



where Ai and A 2 are arbitrary parameters. However, in the particular case K + = 2 the Dirac 
bracket becomes local, since the nonlocality is eliminated due to the contraction of the matrix 
with its inverse matrix. Indeed, one has 



V 



1-A" = (l-A- 1 )^ A fe -^oE A " 

k=u+l k=l 

l-(-l)"A" = (l + A- 1 )^ J2 (-l) fc A fc -^oE(-l) fe A^ 



k=v+l k=l 

and for K + = 2 the Dirac brackets are local 

-2 



A 2 (u n , i ,u md ) = l/2(-l)% n>i [(l + A" 2 ) 

2-n m-2 



fc=3— n k=l s=m—l s=l 
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As an example, we finish this section discussing the explicit form of the second Hamiltonian 
structure of the ID (2, 2)-Toda lattice hierarchy. The Lax operator defining this hierarchy is 
parameterized as follows: 

L 2 ,2 = u ie 29 + %e 9 + a + Pie' 9 + die' 29 . 

Using eqs. (J39j) . (J4(J|) and (J45)) one can derive the corresponding first Hamiltonian structure 

{di, Cj}i = (-iy di(5ij - 5 iij+2 ), 
{ci,Pj}i = -(-i) 3 Pj(8ij-i + S%j), 
{PuPj}i = (-l) j (di5ij-i -djSij+i), 

{ 7i ,7;}i = (-1)^(^-1 -5 iJ+1 ) (46) 
and the second Hamiltonian structure 

{d u d,} 2 = l/2(-l)^(l + A)(5 M _ 2 -5 M+2 ), 

{d i)Pj } 2 = l/2(-iyd iPj ((l - A)(6 itj + 6 i>j+1 ) - (1 + A)^--! + 5 hj+2 )), 
{di,Cj} 2 = (-1) 3 diCjtfij - S id+2 ), 

{d hlj } 2 = l/2(-l) J 'dj 7j -((l - A)(5ij + 6 iJ+a ) - (1 + A)(5 i>j+1 + 5 iJ+2 )), 
{d h Uj} 2 = l/2(-iydiUj(l - A)(6ij - 5ij_ 4 ), 

{Cj, Cj} 2 = (-l) J (M i c/ i (5j ii „2 - UjdiSij+2 + JiPjSij-t + JjPiSij+x), 

{ci,Pj} 2 = ( 1 )' ! (d i" j'^i.j ■ i + djji5ij-2 - CiPj(5i tj + Sij-i)), 

{ci,ijh = -(- i y{ u iPjSi,j-i+ u jPiSi,j+2), 

{PhPjh = (cidjSij-i - <-jd,<h.j . i - l/2piPj(l + A)(8 itj+1 + 

{pulj}2 = (— 1) J (uidjSij-i - Ujdi5ij + 3 - paj(5ij +1 - 1/2(1 - A)(5 iij + ^j_ 2 ))), 

{pi, Uj} 2 = l/2(-iyp i u j (l - A)(5ij - 5 itj+1 + 5 id+2 - 5i )j+3 ), 

{li,lj}2 = (-1) J (uiCjSij-! - UjCi5i :j+ i - 1/27^(1 - A)(<5ij+i + ^--l)), 

{li, Uj} 2 = l/2(-l)- J 7iM j (l - A)(8 iij - 8ij+i + 5i d+2 - Sij-i)), 

{u u u 3 } 2 = l/2(—i) j UiUj(l - A)(5 it j +2 - 5ij-<z)i (47) 

where only nonzero brackets are written down. Here we have introduced the parameter A which 
for the unreduced brackets is equal to zero, A = 0, and for the Dirac reduced brackets with 
the reduction constraint itj = 1 is equal to one, A = 1. In the latter case, algebras (J4lj|) and 
reproduce, respectively, the first and second Hamiltonian structures of the ID generalized 
fermionic Toda lattice hierarchy found in JT] by a heuristic approach. 



7 Bi-Hamiltonian structure of 2D fermionic K )-TL 
hierarchy 

In this section, we construct the bi-Hamiltonian structure of the 2D fermionic (K + , K~)-TL 
hierarchy. This hierarchy is associated with two Lax operators ()24|) belonging to the operator 
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space (J3J). Following [13] we consider the associative algebra on the space of the direct sum of 
two difference operators 

g:=0+ + ©0-_. (48) 

However, in contrast to the case of pure bosonic 2D TL hierarchy, the difference operators in 
the direct sum ()48|) can be of both opposite and equal diagonal Z 2 -parity. It turns out that the 
Poisson brackets can correctly be defined only for the latter case. In what follows we restrict 
ourselves to the case when both operators in g (|4*S|) have the same diagonal parity. 
We denote (x + ,x~) elements of such algebra g = with the product 

, x l ) • {X21 x 2 ) = {x^Xy , XiX 2 ), (49) 
and define the inner product as follows: 

< (x + , x~) >:= str(x + + x~), (50) 
where x + G , x~ G O^. Using this definition we set the Poisson bracket to be 

{A,/ 2 } =< (Qi + ,o-_),[VA,v/ 2 } e >, (51) 

where 

[v/i, v/ 2 }® := ([va + , (v/ 2 + )* ( V +) } , [ V /r, (vj-rrV-)}), 

fk are functionals on g (|4*H)l . and V/fc[(0^ + , Q^-_)] = (V/^",V/^") are their gradients which 
can be found from the definition 



a/ fc [(0+ + ,0"_) + e(50+ + ,50-_)] 



de 



<(50+ + ,50 i ,_),(V/+,V/ fc -)> 



6=0 



< £Q+ + , V/+ > + < so~_, v/ fc - > . 



In order to obtain nontrivial Hamiltonian dynamics, one needs to modify the bracket (|5T|) 
applying the i?-matrix 

[V/i, V/ 2 }® — > [V/ l5 V/ 2 }® = [i2(V/i), V/ 2 }® + [V/ l5 J2(V/ 2 )} e . 

The i?-matrix acts on the space (|48|) in the nontrivial way and mixes up the elements from two 
subalgebras in the direct sum with each other 

R(x + , x') = (x% - x± + 2xZ,xZ - xl + 2x%) (52) 

which is a crucial point of the i?-matrix approach in the two-dimensional case ^3]- This 
i?-matrix allows one to find two compatible Poisson structures and rewrite the Lax-pair repre- 
sentation ()25j) in the Hamiltonian form. 

By construction the i?-matrix ()52)1 satisfies the graded modified Yang-Baxter equation 

R([(x + , x~), (y + , y-)} R ) - [R(x+, x~), R(y + , y~)} = a[(x + , x~), (y+, y')} (53) 
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with a = 1. In order to show this, we simply repeat the arguments of ^3] representing the 
i?-matrix as the difference R — H — ft of two projection operators 



n(x + , x ) = (xi + x_, x^ + x_), n(x + ,x ) = (xt — X_, X + — X 
u(x + ,x~) + n(x + ,ar) = (x + ,x~), n 2 = n, n 2 = n, nn = nn = o. 

Therefore, the i?-matrix (|52j) provides the splitting of the algebra g and solves the modified 
graded Yang-Baxter equation 

The two-dimensional .R-matrix is not graded antisymmetric, its adjoint counterpart W looks 

like 

R\x + , x ) = (x^ — x~^ + 2a;^ , x >0 — x^ + 2x^ ) — H — W , 
where the dual projections are 



ITT/" ry*~ ^~ rtf \ ( K*~ ^~ I ry Hf ^ I rf 

11 \X ,X ) — \X^ Q -f x^ , X >Q -f X 



The direct verification by substitution in shows that the graded antisymmetric part 

also satisfies the graded modified Yang-Baxter equation Therefore, by Theorem of Sec- 

tion 3 there exist two Poisson structures on g (f4*KJ) . 

Using eqs. (}2"UH2"Tj) . (|mHoT?j) . (|5^j) and cyclic permutations inside the supertrace Q we 
obtain the following general form of the first and second Poisson brackets: 

{f,gh = <^(V 9 + ,V< ? -),(V/ + )^ ) > 

+ <Pr(Vg + ,Vg-),(Vf-y {d ^ >, z = l,2, (54) 

where g?v 9 := c?vg+ = dy g -- The Poisson tensors in eq. are found for any values of 

(K + ,K~) for the first Hamiltonian structure 

P 1 + (Vg + ,Vg-) = l(VgI-Vgl)< K+ \(L + K+ y^} 

- ([L+ + ,Vg + } + [L-_,Vg-}U Q , 
P{{Vg\Vg-) = [(VgX-Vg+)< K -\(L-_)< d ^} 

- ([L + K+ ,Vg + } + [L K _,Vg-}) >0 , 

while for the second Hamiltonian structure we constructed the explicit expression of the Poisson 
tensors for even values of (K + , K~) only 

P+(Vs+,V(T) = l/2([(Vg-(L K _y^ + L K _Vg- 

- Vg + {L+ K+ y^ - L+ + V<7+)_, (£++)*«•>} 
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- L+ + ([L+ + ,V<7 + } + [^-,V<T}ko 

- ([L% + ,Vg + } + [L-.,Vg-}Uo 
P 2 (Vg + ,Vg-) = l/2([(Vg + (L + K+ )*^ + L + K+ Vg + 

- Vg-(L K _Y^ - L K _Vg-) + , (L K _)< d ^} 

- L~_ ([L+ +1 Vg + } + [L-_,Vg-}) >0 

- ([L+ + ,Vg + } + [L-_,Vg-}) >0 (L~_)<^ 



The Poisson brackets for the functions u n ^ and v n ^ parameterizing the Lax operators 
can explicitly be derived from (|54|) if one takes into account that 

Vu n ,s = (V< € ,V^) = (e^ + ) a (-l) i ^,0), 
= (V< ? ,V^) = (0,e^--^(-l) 4 ^). 
In such a way one can obtain the following expressions: 

{u n ,i,Um,j}l = {-l) j {5~ K+ +5^ K+ -l)(u n+rn _ K + )i 5ij +n - K + 

_ (_l\(rn+K+)(n+K++l) r . \ 

I -LJ a n+m-K+,j°i,]-m+K+ I , 

iii n \, r_iv'U+ /7-i^ m+ ^)(™ +jfs:++1 )?/ ./). . 

~ u n-m+K- ,i$i,j+n-K+^ + (^ ni A'+ — -0 \ V m-n+K+ ,i$i,j+n-K+ 
I 1 J y m-n+if+ ,j u i,j+m-K- 

{Vn,i,V m ,j}l = (-1)^(1 - C,A- - ^m,A-)(^n+m-A-.i* 



_ (_-\\(m+K-)(n+K-+l) 

for the first Hamiltonian structure and 



-K~,jO 



i,j+m—K~ 



n+m 



J2( 5 U - 5 m , k )((-l) mk u n+m - k ,u k j5, 

k=0 

/ -i \m(n+fc+l) r \ 



i,j+n—k 



— ( — 1)' ? 2 U n,i v m,j[3i,j + $i,j+n-K+ 



m— 1 



~)~ 2 ^ ^ (^Un—m+k,iVk,j+m~k^i,j+n—K+ 



fc=max(0,m-n) 
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(55) 



{Vn,i,V m> j} 2 = (-l)^[f«,iVjft,i-n+X- " (-^^j+m-J?-) 
n+m 

- E(^m,fc - $m,k) ((-l) mk Vn+m-k,iVk,j8 itj - n+k 
k=0 

- (-l) m{n+k+1) v k ,v n+m „ k J lJ+m „ k )] (56) 



for the second Hamiltonian structure; the latter is valid for even values of (K + , K~) only. 

The reduction, according the reduction constraint uq^ = 1, does not require any correction 
terms for the first Hamiltonian structure (|B3j) because {uo t i,u n j} = {uo t i,v n j} = 0. For the 
second Hamiltonian structure the correction terms are 

A 2 (u n>i ,u m>j ) = l/2(-l)V.(l-A x+ - n )(l + A^ + ) 

(A A ' + - A- K+ )-\l + A K+ )(1 - (-l) m A m - K+ )u m , 3 , 

A 2 {u n>i ,v m>j ) = l/2(-l)% n ,(l-A x+ -")(l + A" A ' + ) 

(A K+ - A- K+ )-\l + A K+ )(1 - (-l) m A K - m )v m ,j, 

Mvn,i,v md ) = l/2(-iyv n ,i(l-A n - K ')(l + A- K+ ) 

(A K+ - A- K+ )-\l + A K+ )(1 - (-l) m A K ~~ m )v mij (57) 

and they are nonlocal due to the presence of (A K+ — A _x+ ) _1 in the r.h.s. of eq. (|57|) . 
However, there are unique values of (K + , K~) when nonlocal terms are eliminated. Indeed, for 
K + = K~ = 2 eqs. (JSTj) become local 

A 2 (u nii ,u mtj ) = l/2(-l) J M nii (l + A~ 2 ) 

2-n m-2 



«2 E Afe -c, 2 E Afc )(^ 2 E (-^-^E^ 1 )^ 



k=3—n k=l s=m—l s=l 

A 2 (u nih v m>j ) = l/2(-l)\ ! (l + A- 2 ) 

2-n 2-m 

« 2 E Afe -^ 2 E Afc )(^ E (- 1 ) SAS -^, 2 E(- 1 ) SAS )^' 

fc=3— n k=l s=3—m s=l 

A 2 (v n7i ,v mJ ) = l/2(-l)X i (l + A- 2 ) 

n-2 2-m 

E Afc -^ 2 E Afc )(^ 2 E (-o^-^Et- 1 ) 8 ^ 

k=n~ 1 fe=l s=3— m s=l 

The Hamiltonian structures thus obtained possess the properties (jl2fJTlJ) with do = d L 
d L - . Using them one can rewrite flows (}2"TH30j) for even values of (K + ,K~) in the bi- 
Hamiltonian form 

n± ( U n,i \ _ (( U n,i \ tt± \ _ J ( U n,i \ ri - r - \ 
U »\ (r) -{[ ? >) J , ^s+lll — 1 I (r) j^sh, 

n.i / V n,i / V y n,i 
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with Hamiltonians 



1 1 

Ht = -str(L+ + )l = - (-iNW 



j=— oo 



H; 



i i 00 



For odd values of (fC^K - ) one can reproduce the bosonic flows of (}2?H30|) only. In this case 
eqs. (158)) . due to relation (|17)l . give only fermionic nonzero Hamiltonians using which the 
bosonic flows can be generated via odd first Hamiltonian structure (JHHj) 

n± Z' U n,i \ _ f ( U n,i \ tt± \ 
U 2s\ (r) J ~ i I (r) j^2s+lh- 

n.i / V n.i 



One remark is in order. In Sec. 6, we derived the bi-Hamiltonian structure for the ID 
fermionic (K + ,K~)-TL hierarchy in the i?-matrix approach applying the i?-matrix formal- 
ism developed in Sec. 3. However, the ID fermionic (K + ,K~)-TL hierarchy is obtained as 
a reduction of the 2D fermionic (K + ,K~)-TL hierarchy with the reduction constraint (J3EJ). 
Therefore, this reduction constraint can be carried over into Hamiltonian structures and the bi- 
Hamiltonian structure for the ID fermionic (K + , K~)-TL hierarchy can equivalently be derived 
from that of the 2D hierarchy just by reduction with the corresponding constraint. Actually, 
this reduction amounts to the extraction of subalgebras in the Hamiltonian structures for the 
2D fermionic (K + , K~)-TL hierarchy. Indeed, one can verify that the fields u n ^ (0 < n < K + ) 
and v n> i (0 < n < K~ — 1) form subalgebras for even values of (K + , K~) in both the first ()55)) 
and the second (|56j) Hamiltonian structures and these subalgebras are the first (|39p and the 
second ()40)) Hamiltonian structures, respectively, if one redefines the fields as follows: 

U n ,i = U n ,i, < n < K + 

u K++K -_ n i = -v ni , < n < K~ - 1. 



8 Conclusion 

In this paper, we have generalized the i?-matrix method to the case of ^-graded operators 
with an involution and found that there exist two Poisson bracket structures. The first Poisson 
bracket is defined for both odd and even operators with Z 2 -grading while the second one is 
found for even operators only. It was shown that properties of the Poisson brackets were 
provided by the properties of the generalized graded bracket. We have deduced the operator 
form of the graded modified Yang-Baxter equation and demonstrated that for the class of graded 
antisymmetric i?-matrices it was equivalent to the tensor form of the graded classical Yang- 
Baxter equation. Then we have proposed the Lax-pair representation in terms of the generalized 
graded bracket of the new 2D fermionic (K + , K~)-Tod& lattice hierarchy and demonstrated that 
this hierarchy included all known up to now 2D supersymmetric TL equations as subsystems. 
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Next we have considered the reduction of this hierarchy to the ID space and reproduced the 
ID generalized fermionic TL equations ^Tj. Finally, we have applied the developed R-matrix 
formalism to derive the bi-Hamiltonian srtucture of the ID and 2D fermionic (K + ,K~)-TL 
hierarchies. For even values of (K + ,K~) both even first and second Hamiltonian structures 
were obtained and for this case all the flows of the 2D fermionic (K + , K~)-TL hierarchy can be 
rewritten in a bi-Hamiltonian form. For odd values of (K + , K~) odd first Hamiltonian structure 
was found and for this case only bosonic flows of the 2D fermionic (K + , K~)-TL hierarchy can 
be represented in a Hamiltonian form using fermionic Hamiltonians. 

Thus, the problem of Hamiltonian description of the fermionic flows of the 2D fermionic 
(K + , K~)-TL hierarchy is still open. Other problems yet to be answered are the construction 
of the second Hamiltonian structure (if any) for odd Lax operators and of the Hamiltonian 
structures (if any) for Lax operators L\ and of opposite Z 2 -parities. All these questions 
are a subject for future investigations. 

Acknowledgments. We would like to thank A. P. Isaev, P.P. Kulish, and A. A. Vladimirov 
for useful discussions. This work was partially supported by RFBR-DFG Grant No. 04- 
02-04002, DFG Grant 436 RUS 113/669-2, the NATO Grant PST.GLG.980302, and by the 
Heisenberg-Landau program. 



Here we show that the graded modified Yang-Baxter equation (JTHj) for the case of graded 
antisymmetric operators R is equivalent to the tensor form of the graded classical Yang-Baxter 
equation introduced in the pioneer paper [T3] 



Let Q be a superalgebra with the generators e M (fi = 1, . . . , n + m), the structure constants 
C p and the graded Lie bracket 



where is the Grassmann parity of the generator e M and — 0, if /J, — 1, . . . , n and d^ — 1, 
if // = n + 1, . . . , n + m. 

For the generators of the algebra Q the graded modified Yang-Baxter equation (|T9"j) at a = 1 
takes the form 



Appendix A 



[ru, rn + r 23 ] + [r 13) r 23 ] = 0. 



(A.l) 



{e M ,e^} = e^e u - (-1) 




R({Re^, e 7 }) + R({e^, _Re 7 }) — {Re^, _Re 7 } = {e^, e 7 } 



which can equivalently be rewritten as follows: 



(R^C^R 1 ^ + R^C^R 1 ^ — RPR^C^p) e M — C^e^. 



(A.2) 
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Let us define an invariant supersymmetric non-degenerate bi-linear form associated with 
some representation of Q 

< x, y >:= str(x y) for any x,y G Q 
using which one can introduce the dual basis e M in superalgebra Q 

<e\e v >=5» 

and the supermetric 

^ =< e M , e v >, tT =< e", e v >, jf v = -(-1)V = 

by which one can raise and lower indices as follows: 

e a = Vape?, e a = (-l) d " V a % = e^ a , ^C^yf* = C%. (A.3) 

Note that supermetric t]^ and R£ are even matrices, i.e., for any their nonzero entry one has 
dfj, + d v = 0. 

Now we take the graded tensor product of both sides of eq. (jA.2|) with (g> e 7 and using 
relations (jA.3|) rewrite it in the following form: 

((-l^irtT^C^ + (-l) d >*R xl3 R a »C% a 

-(-!)*> B?#BrCZ f) )ep®e v ®e x = ^V 7 ^ e M ® e„ ® e A , (A.4) 

where 

= (-1)^7/^, Ri = v ai R^. 

Here we use the graded tensor product 

(e„ ® e„)(e A ® e ? ) = (-l)^ dA (e M e A ® e„e s ), (A.5) 
(e„ ® e,)}^ = (e M ), m (e,)^ 2 (-l)* 2 ^ + ^), (A.6) 

where di(dj) means the Grassmann parity of the row (column) of the supermatrix element (e M )^- 
[TH] and one has cfj + cfj = d M for any nonzero (e^)^. 

Assume further that the operator i? is graded antisymmetric, i.e., 

= -(-1)^^ ^ <Rx,y>=-<x,Ry> 
and define two (n + m) x (n + m) matrices 

r = R aP e a ®e P) t = {-l) da r]^e a ® e p , 

where f is antisymmetric, rl 112 = —r^tV: and t is the tensor Casimir element invariant with 
respect to the adjoint action 

[£, z <g> 1 + 1 £3 x] — for any x G 
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Note that r and t are even matrices, i.e., d^ + d i2 + dj 1 + dj 2 = for any rV[*? 7^ or t^ 2 7^ 0. 
Defining triple graded tensor products 

fia =R* f *{e a ®ep® 1), fi 3 = -R a/3 (e a ®l®e /9 ), r a =i{' , ' i (l®e a ®e / ,), 
<i2 = / a (e a ®e^®l), t 13 = / a (e a ®l® e/3 ), t 23 = r] Pa {\ <g> e a <g> e/3 ) 

and using eqs. ()A.2|) and (jA.5j) one can find that eq. (|A.4|) reads 

[ri2, fia + f 23 ] + [fi 3 , f 23 ] = — [ti2, ti 3 ] (A.7) 

that is the tensor form of the graded modified classical Yang-Baxter equation (see e.g. [23]). 
In order to reproduce the graded classical Yang-Baxter equation with the zero in the r.h.s one 
needs to introduce a new matrix 

r = f + t, fi2 = l/2(ri2 - 7*21), £12 = l/2(ri2 +r 2 i) 

for which eq. (jA.7|) takes the form (jA.lj) . Thus, for the case of graded antisymmetric operators R 
we have established the equivalence of equations (J19|) at a = 1 and (jA.lj) which are, respectively, 
the operator form of the graded modified classical Yang-Baxter equation and the tensor form 
of the graded classical Yang-Baxter equation. Note that the former equation is a more general 
one, since it admits solutions which are not graded antisymmetric. 
For completeness we give here the component form of eq. (jA.lj) 

r hi2 r k )dj 2 (d i3 +d j3 ) _ J\iz J* ii / i \d i2 (d i3 +d j3 ) 

k h his \ ' k j 3 jij 2 \ 1 

4. r W3 r iik /-i \di 2 (d il +dj 1 ) _ iiiz i\ k /_i \d j2 (d il +d jl ) 
~ jik'h33\ > hk' jij3\ 1 

1 r hi2 r k i 3 _ 1213 r h k _q 

~ iifc hh k j 3 jij 2 w ' 
where eq. (jA.6j) is used when calculating the triple tensor products. 



Appendix B 

In this Appendix we investigate bi-linear bracket (|21|) and establish conditions on the R- 
matrix and its graded antisymmetric part which are necessary in order the bracket (j2*Tj) be the 
Poisson bracket. Therefore, we need to verify the Jacobi identities for any /, g and h in g 

(-l) d ^{h,{f,g} 2 } 2 + c. p. = -l/4(-l)^<[0,V({/^} 2 )}- 

R^iy h)*( d9+df ^o* ( - df+d9+dh ' > + o* (d9+d/ )(v/i)* (ds+df ^ 

- R(v({f,g} 2 )0< d * +d f)+OV({f,g} 2 ))- 
P*(d g +d } )^ vh y(d g +d } )y > + c p = 0> ( B1 ) 

where 

V({/, g} 2 ) = -1 /4 [vp R((yfy {d3) 0* {d f +ds) + 0* {ds) (V/)* (ds) ) 
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+Vg R}((Vf)* {da) 0* {df+da) - 0* (ds) (V/)* (ds) ) 
+R^{OVg - VgO*^) (V/)* {d9) ) 
_(_iy f d a ( R (yf®*(d f ) + V/) (Vg)* {df) 

+V/ R((\/g)* {df) 0* {df+d9) + 0* {d f ) (yg)* {df) ) 
+Vf R j ((Vg)* {df) 0* {d f +d3) - 0* (c!f) (V#)* {(J/) ) 

+R^{VfO* (df) - OVf) {Vg)* (d f : 



(B.2) 

Inserting (|B.2|) into ()B.1|) after tedious but straightforward calculations we get for the Jacobi 
identities 

(-l) d/(dh+ " 9) < [0,/i}([ J R t F_,i?t G _} + [ RF+ ,RG+} - R([F + ,G + } R )^ > +c.p. = 0, 

where we introduce the notation 

F± = Q*( d h) (y fy( d h) ± (y fy(dh)Q*(d h +d f ) 

G± = ®*( dh+d f"> (\7 gy( dh+d f1 ± (ygy( d h+df)®*(d h +d f +dg) 

for brevity. For any linear map R and its graded antisymmetric part A = 1/2(R — R^) one has 
the identity 

(_iyAd h +d B ) < [o,/ i }[i?tF_, J RtG_} > +c.p. 

= (-i)d f (d h +d a ) < [tyh}U/3([AF-,AG-} - A([F + ,G+} A ) 
- ( [RF_ , RG^} - R([F_, G_} R j) > +c.p. (B.3) 



Now using ()B.3|) and the following identity 

^d f (d h +d g ) < [0; ^ G _| + 3[F+; > +c p = 0> 

which can directly be verified, we finally rewrite the Jacobi identities with an arbitrary param- 
eter a as follows: 

(-1)^*+^ < [0,h}([RF+,RG+} - R([F + ,G + } R ) +a[F + ,G + } 
- ([RF-, RG-} — R([F~, G-.}r) + a[F-, G-} 
+ 4/3([AF_, AGL} - A([F_,G_} A ) + G_}))> +c.p. = 0. 



Now it is obvious that the Jacobi identities are satisfied if R and its graded antisymmetric part 
A obey the graded modified Yang-Baxter equation ()19|1 with the same a. 
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